A self-similar model for shear flows in dense granular materials 
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We propose a model to describe the quasistatic shearing of dry granular materials, 
which notably captures the differences in velocity profiles recently observed in 2 and 3-D 
Couette flow experiments. In our scheme, the steady-state flow is due to the intermittent 
motion of particle clusters moving together with the wall. The motion of a cluster is 
associated with the transient formation of a fracture inside the sheared pack. The model 
is based on the existence of a persistence length for the fractures, which imposes a self- 
similar structure on the clusters. Through a probabilistic approach, we can evaluate the 
rate of appearance of a cluster of a given size and obtain a prediction for the average 
velocity profiles. We also predict the existence of large stress fluctuations at the moving 
wall, which characteristics are in good agreement with experimental data. 



I. INTRODUCTION 

Dry granular materials exhibit a large range of dynamical behaviours (for a general review see ffl). 
An assembly of non-cohesive macroscopic particles confined in a box can be considered as a solid, but 
will start flowing when submitted to a large enough shear stress. However, the underlying mechanism 
can be very different, depending on the shearing rate imposed and the resulting density of the material. 
Upon rapid shearing or shaking, a fluidized state is reached in which the stress is transfered from the 
boundaries to the bulk through binary collisions. This dynamics is well understood under the scope of 
the so-called kinetic theory In many practical situations however, gravity maintains a high density 
in the material so that the stress is mainly transfered by rubbing friction between particles in persistent 
contact. In this regime, all relative motion between particles tend to be confined to very narrow regions of 
the material, making a continuous approach inadequate. Within these so-called shear band, each particle 
is submitted to large forces which rapidly fluctuate as the packing environment and the associated force 
network evolve. The flow is driven by collective and jerky moves of large sets of particles, hold together 
by transient force chains. Although this regime is relevant to soil mechanics and geology (earthquakes 
and pyroclastic flows for example J^j) and important in industrial applications (granular transport in 
hoppers and pipes), a satisfactory description is to date lacking for such fully developed dense flows 

A strong experimental effort has however been seen recently to probe the grain-scale dynamics within 
these shear bands, which now provides a good test for physical models: two series of experiments have 
been performed independently on 2-D and 3-D Couette cells. A collection of particles is confined between 
a fixed outer cylinder and a rotating inner one. In the 2-D experiment, the system consists of one 
monolayer of disc- like particles squeezed between 2 horizontal plates. In this geometry, a shear band is 
always present at the vicinity of the rotating inner cylinder. By different techniques, the authors were 
able to precisely measure the decay of the average particle velocity from the inner moving wall to the 
immobile region towards the outer cylinder |^||. Surprisingly enough, the two results were found to be 
significantly different: in 2-D the velocity V decayed with the distance r to the inner wall according to an 
exponential law, whereas in 3-D, the velocity profile was well fitted by a Gaussian centered on the surface 
of the wall. These results showed to be robust for polydisperse or highly irregularly shaped particles. For 
round and monodisperse grains however, a layering effect and a strong slippage between adjacent layers 
were observed which result in more complicated profiles. 

In this note we will focus on the Gaussian/exponential behaviours as the spatial dimension changes. 
Because a Gaussian form cannot be derived from a local differential equation (for example see fit]), the 3-D 
velocity profile suggests that a purely local model would fail to describe it. Here we introduce non-local 
effects by postulating a correlation length in the particles displacements which depends on their position 
within the shear band. This approach comes to modeling the flow as a succession of transient fractures 
allowing the coherent motion of clusters of various sizes. This hypothesis alone allows us to produce an 
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average velocity profile consistent with what was observed in both 2- and 3-D. Finally, we will show that 
the probability distribution and the force spectrum derived from this model are in good agreement with 
experimental results found in the literature. 

II. THE MODEL 

We consider a plane shear geometry: a wall is moving at a constant speed Vq along a half-space of 
granular material. The flow velocity is assumed to vanish far from the moving wall (see figure |l|) . We 
impose a no slip condition at the wall so that V(0) = Vq (this boundary condition mimics a classical 
experimental realization where the first layer of particles is glued onto the moving wall). In our scheme, 
the motion of each particle can be decomposed into successive steps of typical length d g occurring at 
velocity v = Vq. Between two steps, the particle remains immobile. One grain can thus either be at 
rest or move together with the wallf]. The interface between moving and static particles defines a failure 
surface for the pack. We define P(r) as the stationary probability of a failure surface to be present, per 
unit of time and length (in units of particle diameter) along the y-axis, at a distance r from the moving 
wall (as in ref 0). Although a reduction in granular density by roughly 10% - the so-called Reynolds 
dilatancy - is known to accompany shear flow, we will first consider a uniform density. Both two and 
three dimensional geometries will be addressed; D will denote in the following the space dimension. A 
'fracture' will either be a line in 2-D, or a plane in 3-D. 
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FIG. 1. Self-similar construction describing the shear flow. The wall moves at a constant velocity Vq as well as 
the first layer of particles. At a distance r from the wall, the flow is due to the sudden motion of clusters of size 
r (or bigger) in all directions. The lines represent the possible location of fractures that appear in the deforming 
pack during the shearing process. 

We first consider the dynamics of the first freely moving layer in the region d g < r < 2d g . As the 
wall moves by a distance d g in a time tq = d g /Vo, a particle located in this layer is cither dragged 
by the wall producing a failure further away, or stays immobile so that a crack develops at exactly 
?' = d g . Assuming a Coulomb-type slip condition, the probability of the latter depends on the ratio 
of the normal and shear stress to which is submitted the particle during the time period tq. Based on 
this mechanical consideration, and by analogy with thermally activated processes, Pouliquen & al 
have proposed an expression for the probability po of slippage between 2 particles from 2 consecutive 
layers in a granular material slowly flowing down a pipe. In their picture, the spatial stress fluctuation 
induced by the randomness of the packing structure plays the role of a temperature, allowing them to 
apply classical rate processes theory [G3. In the quasi-static regime we address here, the average shear 
and normal stress are uniform in the material since each layer is constantly in mechanical equilibrium. 
In the limit of uniform density, it is therefore natural to assume a uniform value for pq. Although, the 
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derivation of the slippage rate not only requires the knowledge of the static probability po. A typical 
fluctuation time-scale, related to the flow itself also needs to be introduced. This time-scale corresponds 
to the relaxation time of the stress network around the chosen particle. For a particle in the first layer, 
the only relevant time-scale is tq which is the time needed for a shear strain of 1 to be established. We 
can now write down the rate of slippage at the wall P(d g ) as : 

, n ■ dr Vq ■ dt dr 

P{dg)-dt- _=^.p .-L_ (1) 

a g a g a g 

A is a constant introduced for normalization. It should be noted that we implicitly suppose a clear 
separation between the successive hopping events or equivalently that po << 1 . 

The calculation of P(r) away from the wall is based on a self-similar argument. We postulate that, 
for a yielding event to occur at a distance r from the wall, the crack must extend radially over a length 
of the order r. Hence the shearing of the material at a distance r from the wall requires the coherent 
motion of a solid cluster of size r in all directions. This cluster is dragged by the wall over a distance 
d g as one solid object, in the same way as was the single particle from the first layer (see figure Q). 
This self-similar description of the cluster size is introduced to account for the perturbation of the force 
network in the vicinity of the wall. Chain forces, which are responsible for the rigidity of the dense 
pack, are screened by the proximity of a solid boundary. Independent motion of particles are thus only 
allowed in the vicinity of the wall whereas the bulk behaves as a solid body. We note that an identical 
scaling have been successfully postulated in (at least) 2 very different contexts to describe the effect of a 
wall on dynamical structures. To estimate the Prandtl mixing layer in turbulent flows near a wall, the 
characteristic length of the vortices are assumed to increase like the distance to the boundary layer. The 
derivation of polymer dynamics near a wall in the semi-dilute regime also assumes that the 'blob' size 
(the dynamical coherence length of the monomers) grows like the distance to the wall |llj]. As the width 
of the shear band in granular flows is always of the order of 10 bead diameters, regardless of the bead 
size, the bead diameter appears to be the only relevant length-scale in this problem. A linear increase of 
the coherence length with the distance to the wall is therefore the only reasonable choice^. For simplicity 
we will take exactly r as the coherence length in the following. 

Given this general picture, the scaling analysis follows: we define N g (r) as the number of particles 
involved in a crack developing at a distance r : 

/ \ D-l 

The probability for this set of particles to become simultaneously unstable is {po) N9 , whereas the 
time-scale r(r) of the stress fluctuations experienced by the cluster of size r reads: 
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The scaling for r(r) was made consistent with the criterion used for the first layer: r(r) corresponds 
to a shear strain of 1 for the cluster of size r. Two different time-scales are thus involved in the dynamics 
of these blocks: a rapid one, tq associated with the release of the stress when a yielding occurs; a slower 
one, r(r) which characterizes the stress fluctuations experienced by the particles at a distance r from the 
wall. 

Finally, a particle located at y = r has N g possible locations along the crack and the expression for 
P(r) eventually reads: 

n , n , dr , „. , .»T Vndt dr . . 

P(r)dt ■ — = A ■ N g ■ (po)"» ■ • — (2) 

dg r dg 

When yielding occurs at a distance r from the wall, every particle below y = r moves over a distance 
d g whereas no motion occurs above y = r. Such an event represents a discontinuity in the instantaneous 



2 A different mapping of the coherence length Z(r) may have to be used for a different geometry. For example, 
the stationary dense granular flow along an inclined plane seems to exhibit a characteristic length H depending 
on the inclination and roughness of the inclined plane jl2]j]l|. We postpone the discussion of the coherence length 
in this case for a further study. 
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velocity profile (the velocity is Vb below the crack and above during a time to). The constitutive 
differential equation for the stationary mean velocity profile V(r) thus reads : 



dV 

= -P<r) (3) 



with the boundary conditions V(0) — Vq and V(oo) = 0.. 
For D = 2, we obtain: 

0V 2 D A r/d 

— = -d- g Yoipo) 

which by integration gives an exponential velocity profile V2d{t)'- 



For D — 3, it gives: 



Vid =Vo ■ e Ad s with A = -!— and A = — 

ln(Po) ^ 



dV 3D _ Ar {r/d 



so that the velocity profile V^oir) is Gaussian: 
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V 3D =Vo-e 2{ " d ^ 2 with cr 2 = - - . , and A = — 

21n(p ) ct 2 

The function P(r) then reads respectively in two and three dimensions: 

P2D(r,t) = rrjj-e ^ and P 3D (r,t) = -j -^e 

Although a few hypothesis have been used to derive these results, one should note that the major 
assumption is that the persistence length of the fractures increases linearly with r. This argument alone 
produces the main behaviour of P(r) as a function of pq. The rescaling of the stress fluctuation time-scale 
r(r) only controls the pre-factor of the exponential and Gaussian velocity profiles. A correction could 
also arise from a spatial dependence of pq. For instance, the shear induced density profile neglected in the 
present approach might affect the uniformity of the average shear and normal stress per particle within 
the shear band. Considering that po is a decreasing function of the density, this correction may increase 
the velocity decay and narrow the shear band. This may explain deviations from the Gaussian and 
exponential profiles observed in experiments with rounded and monodisperse particles where the largest 
density gradients are present. However, an exact solution of (^) for a general po is highly non-trivial and 
should be investigated in a further work. 

III. FORCES DISTRIBUTION AND POWER SPECTRUM 

In addition to the velocity profiles, the knowledge of the function P(r) gives appropriate information 
for computing other properties of the flow such as the probability distribution function (PDF) of the 
forces acting on the wall as well as the power spectrum of this signal. These two quantities have been 
experimentally studied in both two and three dimensional granular shear flows |^,|l5fl. In the latter 
experiment, the normal force on the bottom of a shear cell was monitored using a force transducer which 
could accommodate a large number (from 4 to 100) of particles. The resulting measurement corresponded 
to the integrated value of many individual contact forces. The measured force signal F(t) appeared as a 
series of narrow peaks of various heights. To compare these results with the present model, we suppose 
that each peak is the signature of a failure in the structure of the sheared pack. Assuming that the 
slippage of one grain releases a given force, the motion of a cluster of size r produces a peak of intensity 
F(r) proportional to N g , or equivalently: 



F(r) oc r 



These peaks last for a time tq and are uncorrelated in time. ^From here, we can derive the form of the 
probability distribution function (PDF) p(F) by identifying r D ~ x with F in P(r): 



P2d{F) = P3d{F) = (4) 

where Fo is the mean force of the distribution. Surprisingly enough, these PDF's are the same for 2- 
and 3-D. Both agree with experimental observations at large forces ||. The behaviour at low forces is 
not probed by experiments although it may have more complicated features. 

One may note that eq. (^) is formally identical to what has been observed for large forces on individ- 
ual particles in a compressed static array []l6|| . But this static distribution would narrow upon spatial 
averaging. By contrast, the width of the force distribution in the present description is insensitive to the 
integration area^ in agreement with experimental observations. In the continuously sheared regime, the 
wide force distribution arises from the coherence in the force release induced by the sudden motion of 
large clusters, and not from the purely static force distribution on individual particles. 

The power spectrum of the force fluctuations can also be obtained in this description as follows: the 
power signal is decomposed into a stochastic succession of peaks of height F 2 (r) and width r . The period 
r'(r) between two successive pulses of amplitude F 2 (r) is simply the inverse of the rate P(r). For a given 
r, and therefore a given F(r), the signal is the so-called telegraph noise, for which the power spectrum 
Q r {uj) reads jL7|: 

F 2 (r) r'(r)r 1/T(r) 

Wr\W) — — 



(r'(r)+T ) 2 ^ + (1/T(r)) 2 

with l/T(r) = l/r'(r) + l/r 

Therefore, since failures at different distances from the wall are uncorrelated, the power spectrum Q{lo) 
for the overall signal is given by the sum: 

QH = -J- / Q r (io)dr (5) 
dg Jo 

Q(lo) exhibits qualitatively the same behaviour for 2 and 3 dimensions, dominated by Lorentzian-likc 
functions (see figure (Q) ). For frequencies much larger than lu — l/r , the power spectrum behaves 
like ujo/uj 2 as observed in experiments |l5[ ]; this is a direct consequence of the stochastic description of 
the process. On the other hand, the expression (^|) displays a well-defined non-zero limit as ui — > 0. By 
comparison, the power spectra measured in ]l5[ ] exhibit a continuous, although small increase as oj — > 0. 



3 The finite size should merely show up as a cut of at high forces. 
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FIG. 2. Power spectrum of the forces in the 2D geometry. The curve has a fiat part at low frequencies and 
a 1/tu 2 behaviour for large frequencies. The cross-over between these two regimes occurs around the typical 
frequency loo for which the wall has moved by one grain diameter. 

IV. CONCLUSION-ACKNOWLEDGEMENTS 

We have proposed a stochastic description of slowly sheared granular material that can capture the main 
characteristics of the flow: the average velocity profiles and force fluctuations in 2 and 3 spatial dimensions. 
We want to underline the relatively good agreement with experimental data with regards to the small 
number of ingredients introduced. In particular, we ignored the existence of granular density variations 
along the direction normal to the flow which may induce significant corrections. The crucial assumption 
of this model lies in the self-similar structure of moving clusters that rapidly form and disappear as the 
material flows. This assumption has been qualitatively suggested by the mere observation of 3-D Couette 
flows watched from below through a transparent bottom plate where the jerky and multi-scale dynamics 
is clearly visible. We hope that this tentative description will motivate further quantitative investigations 
to probe the nature of space and time correlations in the motion of neighboring particles in such systems. 

Here we have focused on recent experimental data in order to test this self-similar model but we think 
that this simple description of dense granular dynamics as a super-imposition of coherent moves can 
be extended to many more geometries and processes than just plane shear. The issue is to correctly 
prescribe the coherence length mapping i(r), and describe the probability P(r) for each realization. In 
particular, we anticipate that the existence of multiple relaxation times recently demonstrated in granular 
compaction experiments ]l8| ] could be associated with rearrangements over different length-scales during 
the relaxation of the packing structure. To that extent, this model may offer a clue to understand the 
numerous features th at g ranular material share with glassy liquids below the glass transition (aging, slow 
relaxation, jamming [[L9|). 
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